
% run this output of this program as a script:

% GraphBLAS matrix A: nrows: 3 ncols 5 entries: 10
A = sparse (3,5) ;
    A (1,2) =  (                   0) + (  0.8253972948485638)*1i ;
    A (1,3) =  (  0.8158771738902956) + (  0.6199473944886393)*1i ;
    A (1,4) =  (   0.230335927869842) + (                   0)*1i ;
    A (1,5) =  (  0.6834455780600209) + (                   0)*1i ;
    A (2,1) =  (                   0) + (  0.8145545554415048)*1i ;
    A (2,2) =  (                   0) + (  0.8106511441383925)*1i ;
    A (2,3) =  (                   0) + (  0.6296551667948371)*1i ;
    A (2,4) =  (  0.6749993971837149) + (  0.1028046574093258)*1i ;
    A (2,5) =  (  0.1781985185630041) + (                   0)*1i ;
    A (3,1) =  (  0.6668132827629559) + (                   0)*1i ;
A

% GraphBLAS matrix B: nrows: 5 ncols 4 entries: 13
B = sparse (5,4) ;
    B (1,1) =  ( 0.03470661760412896) + (                   0)*1i ;
    B (1,2) =  (  0.2320846810101108) + (                   0)*1i ;
    B (1,4) =  (  0.4037695773165162) + (                   0)*1i ;
    B (2,1) =  (                   0) + (    0.70501889367129)*1i ;
    B (2,3) =  (  0.2114568817022959) + (                   0)*1i ;
    B (2,4) =  (                   0) + (  0.4423346124115619)*1i ;
    B (3,1) =  (   0.929676618909714) + (                   0)*1i ;
    B (3,2) =  ( 0.03154825424997802) + (                   0)*1i ;
    B (3,3) =  (                   0) + (  0.1598994553484671)*1i ;
    B (3,4) =  (  0.9813322475732387) + (  0.9048904118202962)*1i ;
    B (4,1) =  (                   0) + (  0.9376008381429164)*1i ;
    B (4,4) =  (                   0) + (  0.8269171297987302)*1i ;
    B (5,4) =  (  0.1371733231185191) + (                   0)*1i ;
B

% GraphBLAS matrix C: nrows: 3 ncols 4 entries: 11
C = sparse (3,4) ;
    C (1,1) =  (  0.1765812448145327) + (  0.7923137566352751)*1i ;
    C (1,2) =  ( 0.02573950051864457) + ( 0.01955825802293901)*1i ;
    C (1,3) =  (-0.09912925072343469) + (  0.3049942538704926)*1i ;
    C (1,4) =  (-0.03168916369751931) + (   1.537122526250736)*1i ;
    C (2,1) =  ( -0.6679141057457941) + (    1.24652612056368)*1i ;
    C (2,2) =  (                   0) + (  0.2089101554568301)*1i ;
    C (2,3) =  ( -0.1006815182278427) + (  0.1714177630879029)*1i ;
    C (2,4) =  ( -0.9889148320979797) + (    1.50496183271377)*1i ;
    C (3,1) =  ( 0.02314283361820783) + (                   0)*1i ;
    C (3,2) =  (  0.1547571480233455) + (                   0)*1i ;
    C (3,4) =  (  0.2692389173302374) + (                   0)*1i ;
C
E = A*B
err = norm (C-E,1)
assert (err < 1e-12)
